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Abstract 

A study of signals originating near the lithium- diffused n+ contact of 
p-type point contact (PPC) high purity germanium detectors (HPGe) is pre- 
sented. The transition region between the active germanium and the fully 
dead layer of the n+ contact is examined. Energy depositions in this transi- 
tion region are shown to result in partial charge collection. This provides a 
mechanism for events with a well defined energy to contribute to the contin- 
uum of the energy spectrum at lower energies. A novel technique to quantify 
the contribution from this source of background is introduced. Experiments 
that operate germanium detectors with a very low energy threshold may 
benefit from the methods presented herein. 

Keywords: germanium detectors, dead layer, transition layer 



1. Introduction 

P-type point contact (PPC) high purity germanium (HPGe) detectors 
[H, 0] are of interest for use in astroparticle physics. The small area of the 
point contact (see Fig. [I] for a diagram) results in a sharply peaked weighting 
potential near the point contact as well as a small readout capacitance. The 
localized weighting potential of PPC detectors results in distinct current 
pulses from individual interaction charge clouds, and enables pulse-shape 
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analysis to be used to distinguish between single-site and multiple-site events. 
Backgrounds from gamma rays scattering multiple times in a detector can 
thus be identified [2L y, |4J resulting in this technology being adopted by 
the Majorana [lITa @, II and GERDA ji[ experiments to search for the 
neutrinoless double-beta decay of 76 Ge. 

The small readout capacitance allows PPC detectors to be operated with 
low levels of electronic noise and a correspondingly low energy threshold. 
This makes them attractive candidates for light dark matter searches, such 
as the CoGeNT experiment 0, [hJ which has demonstrated the operation of a 
PPC detector with a threshold of ~ 400 eV. The low energy threshold of these 
detectors has also made them candidates for detecting coherent neutrino- 
nucleus scattering 0, EH, 12 1. Searches for neutrinoless double-beta decay 



using Ge also benefit from the low energy threshold of PPC detectors; x- 
rays at 1.3 keV and 10.4 keV following the electron-capture decay of 68 Ge 
can be used to tag the subsequent (3+ decay of 68 Ga which has a Q- value 
of 2.9 MeV, above the region of interest for neutrinoless double-beta decay. 
Understanding the behavior of these detectors at low energies is critical for 
all these types of experiments. 



A recent publication by the CoGeNT collaboration [lOj presented re- 
sults from a PPC detector used to examine x-rays with energies of 10.4 keV 
from the decays of internal 71 Ge (produced by thermal neutron activation of 
the detector using an americium-beryllium source). X-rays that deposited 
energy near the thick, lithium-diffused, n+ outer contact of their detector 
were shown to result in signals with incomplete charge collection at ener- 
gies between 2.5 keV and 7.5 keV (Fig. 1 in Ref. [13]). Their observation is 
consistent with the existence of a transition layer, between the bulk active 
germanium and the contact dead layer, where a weak electric field results in 
slow and incomplete charge collection. Signals originating in this transition 
layer give rise to characteristically slow charge pulses and contribute to the 
continuum at lower energies. Understanding the properties of energy deposi- 
tions in the transition region is important for determining the active volume 
of a detector as well as for understanding how the incomplete charge collec- 
tion of these events contributes to the continuum in the energy spectrum. 

The existence of the transition layer between the contact and the active 
germanium has been established in the literature, although a comprehensive 
picture is lacking and past research has been somewhat disconnected. Pulses 
resulting from anomalous charge collection were identified in the early devel- 
opment of lithium-drifted germanium Ge(Li) detectors. Soon after the initial 
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characterization [13, [l4| of Ge(Li) detectors in 1963 by Ewan and Tavendale, 
Alexander et al. [15[ observed a class of events with abnormally long collec- 
tion times, and attributed these to interactions occurring between the fully 
compensated (active) part of their detector and the thick (dead) n+ contact, 
where the electric field is too weak to collect charge carriers efficiently. They 
warned that any electron interactions near the n+ contact, including inter- 
actions from electrons produced when high energy gamma rays interact in 
the surroundings of the detector, would contribute to the continuum of the 
energy spectrum. 



In 1967, Tamm et al. [16] identified diffusion as the principal mechanism 



for charge carriers created in the transition layer to enter the region where 
the electric field is strong enough for effective charge collection. The diffu- 
sion coefficient, the depth of the un-depleted region, and the charge carrier 
lifetime were used to estimate the diffusion time to be of order microseconds, 
giving rise to slow pulses and giving enough time for charge loss due to re- 
combination. Tamm et al. showed that the background from slow pulses 
increased at lower energies in their energy spectra. Detailed experimental 
studies by Strauss and Larsen 17J in the same year and by Sakai 18( in 1971 
confirmed, using a variety of detector geometries and collimated sources, that 
the slow pulses originated near the contacts. 

More recent studies of germanium detectors, focused on efficiency mea- 
surements and the detailed modeling of energy spectra, have also examined 
the properties of events near the lithium-diffused contacts. In 1989, Clouvas 



et al. [19] noted that the efficiency measurements for their detector yielded 



a dead layer thickness of 2.5 mm that did not agree with the value of 0.5 mm 
provided by the manufacturer. They additionally found that in order to 
model the energy spectrum from their p-type coaxial detector at low ener- 
gies, they needed to include a transition layer where partial charge collection 
occurred. By considering the diffusion of the lithium in the n+ contact, they 
derived a model for the size of the transition region that enabled them to 
successfully simulate their energy spectrum at lower energies. 

In 2003, Rodenas et al. [20] showed that the lithium-diffused layer changes 
both the active volume of the detector, and the efficiency for detecting low 
energy radiation due to the increased attenuation from the dead layer. They 
pointed out that the thickness of this layer is difficult to measure because 
of the partially active transition region and that the values provided by the 
manufacturers are often wrong. In 2005, Maleka et al. 21] used Monte Carlo 



simulations to show that small changes in the thickness of the lithium-diffused 



4 



layer can result in significant changes in the efficiency for low energy gamma 



rays. Dryak et al. [22] argued that neglecting the transition layer results in 
underestimates of the total efficiency of a detector in their work from 2010. 
They also showed evidence, using a coincidence measurement, that signals 
from the transition layer were delayed in time compared to normal signals, 
which is consistent with diffusion being the dominant transport mechanism 
for charge carriers in the transition layer. The Majorana collaboration 
is performing detailed measurements using coincidences between betas and 
gammas and between gammas to study the precise mechanism of charge 



carrier transport in this transition layer using the detector described in [23 



and another similar detector produced by CANBERRA. Results of these 
studies will be presented in a subsequent publication. 



Studies between 2007 and 2011 by Huy et al. |2j and by Huy |25|, [26 
showed that the depth of the lithium-diffused layer increased as their HPGe 
detector aged, due to diffusion of the lithium atoms over time. They noted 
that this effect is particularly pronounced for periods where the detector was 
kept at room temperature. Using these data, they showed that the counts 
in the lower part of their energy spectrum increased with the depth of the 



lithium-diffused layer. In 2010, Luis et al. [27J observed significant discrep- 
ancies between the dead layer they measured on a Broad Energy Germanium 
(BEGe) detector and that quoted by the manufacturer (CANBERRA (28|); 
they attributed the discrepancy to an increase in dead layer thickness during 
the time between their measurement and the one performed by CANBERRA. 

In summary, there is substantial evidence in the literature to support 
the claim that high energy events can contaminate the lower part of the en- 
ergy spectrum continuum in germanium detectors that have a thick, lithium- 
diffused, contact. These events occur in a transition region between the fully 
active germanium and the dead outer contact. Pulse shapes from these events 
are characteristically slower than those originating in the active region, de- 
layed in time and have a degraded pulse height. These observations are 
consistent with the hypotheses that the dominant transport mechanism for 
charge carriers created in the transition region is diffusion into the bulk re- 
gion, followed by drift in the electric field. Furthermore, the thickness of 
the lithium-diffused region can change in time due to continuous diffusion of 
lithium. 

This paper presents the results of a study that characterized events orig- 
inating near the lithium-diffused contact of a custom-made PPC detector 
illuminated with a collimated source of low energy gamma rays. The ex- 
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istence of a transition region between the fully active germanium and fully 
dead layer of the n+ lithium-diffused contact is verified. Energy depositions 
in this region are observed to result in pulses with slow rise times and de- 
graded charge collection. A novel method for determining the depth of the 
transition layer and the depth at which energy depositions become fully de- 
tectable is introduced. The method is also applied to the detector described 
in 23] and it is shown that the size of the transition and dead layers is depen- 



dent on the process that was used to fabricate the lithium-diffused contact. 
2. Description of the LBNL PPC detector and data acquisition 



A custom PPC detector [29(, illustrated in Fig. [TJ was fabricated at 



Lawrence Berkeley National Laboratory from a high purity germanium crys- 



tal purchased from ORTEC [30]. The crystal has a radius of 31mm and a 
height of 50 mm. A lithium-diffused layer covers the outside of the entire 
crystal except for the surface with the point contact, where it only covers 
the outer 5mm (the "lithium wraparound"). The lithium-diffused surface 
(except for the wraparound) was metalized by vapor deposition of aluminum 
to create a low resistance high voltage outer contact. A hemispherical dimple 
with a diameter of 1.0 ± 0.5 mm was machined into the crystal, sputter coated 



with amorphous silicon [31] , and metalized by vapor deposition of aluminum 
to define the point contact. The detector was operated with an un-passivated 
surface between the point contact and the lithium-diffused wraparound. 

The crystal was installed in a custom copper mount prototype designed 
for the MAJORANA DEMONSTRATOR experiment 0, @, @, @). The mount 
was fixed to a copper cold plate and surrounded by a 1.5 mm thick copper 
infrared shield to keep the detector cold. The detector assembly was housed 
in an aluminum cryostat with a wall thickness of 2.2 mm and cooled with 
liquid nitrogen. Charge pulses from the point contact were processed with a 
cold FET mounted on a low mass board designed for low-noise performance 



and low radioactivity levels [32] . The amplification loop was closed outside 
the cryostat in a custom charge sensitive resistive-feedback preamplifier with 
a 16 MHz bandwidth. Charge pulses processed by the preamplifier were 
digitized using a Struck SIS3302 16-bit 100 MHz digitizing card. 

The detector was scanned with 59.5 keV gamma rays from a 1 mCi 241 Am 
calibration source collimated using tungsten collimators (10 mm thick, hole 
diameters between 0.5 mm and 2 mm). Two scans were performed; one in 
the xy plane (illuminating the point contact surface of the detector) and one 
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Figure 1: Schematic diagram of the LBNL PPC crystal, showing the point contact (A), 
the 5 mm lithium-diffused "wraparound" (B), the lithium-diffused layer (C), and the alu- 
minized low resistance contact for the high- voltage (D). The coordinate system referred to 
in the text is also shown, with the y axis going into the page and the origin at the point 
contact. The aspect ratio of the detector is correct, although the size of the point contact 
(1 mm diameter), and the thickness of the lithium and metalization layers are not to scale. 
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in the yz plane (illuminating the side of the detector) . Readers are referred 
to Fig. [TJfor the definition of the coordinate system used in this paper. Two 
different setups (collimators, source positioning system) were used for the 
two scans. Consequently, the beam spot sizes on the crystal were different 
in the two scans. The diameter of the beam was determined from geometry 
to be 3 mm ± 1mm in the xy scan and 2.0 mm ± 1,0 mm in the yz scan. In 
each scan, data were taken at various positions of the collimated source for 
a fixed duration of 4 minutes. 



3. Determination of the full charge collection depth 

The thickness of the lithium-diffused layer of a detector is difficult to know 
accurately from the manufacturing process. The n+ contact on the LBNL 
PPC detector was fabricated by first evaporating lithium onto the surface of 
the crystal (masking the areas where no lithium was desired). The crystal 
was then heated up to approximately 280°C in order to increase the diffusion 
of lithium atoms into the crystal. Variables such as oxygen concentration 
at the surface, and the temperature profile of the crystal as a function of 
time during the diffusion process result in large uncertainties when trying to 
determine the lithium thickness from the manufacturing parameters. 

A technique is presented here to characterize the thickness of the region 
where the lithium diffusion affects energy collection, usually referred to as 



m 



the "dead layer" . The method is similar to that presented in [33j and [34 
which a collimated beam of gamma rays was incident on the detector surface 
at different angles. 

It will be apparent from this work that the term "dead layer" must be 
defined appropriately, and that there are three principal regions that must be 
distinguished. The outermost region of the contact, which is conductive and 
has no electric field, may be considered "dead". There is then a transition 
region where energy depositions result in incomplete charge collection. The 
transition layer depth (TLD) is defined as the depth at which this region 
begins and corresponds to the boundary with the dead layer. Finally, the 
full charge collection depth (FCCD) is defined as the depth at which the 
truly active part of the germanium detector begins and corresponds to the 
boundary of the transition region and the region where full charge collection 
occurs. Higher energy depositions can result in charge carrier clouds that 
extend across the different regions. In this case, the different parts of the 
charge cloud will be collected differently, depending on the region in which 
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they were created. The 59.5 keV gamma rays used in this work result in 
electrons with a range of less than 20 /jm in germanium, and the associated 
charge carrier clouds can effectively be taken as point-like. 

Fig. [2] shows a bottom- view of the setup that was used to scan the side of 
the detector in this work. The beam from the source was collimated in the 
x-direction. The copper infrared shield and the aluminum cryostat are also 
shown. It can be seen that, as the y position of the source is varied, the path 
length of the beam through the lithium-diffused layer, xu(y), changes. Given 
the attenuation coefficients [35j for copper, /xc u = 1.63cm 2 /g, aluminum, 
/iAi = 0.281 cm 2 /g, and germanium, fiQ e = 2.07cm 2 /g, the beam attenuation 
from the various materials can be determined. The attenuation coefficient of 
the lithium-diffused layer is taken as equal to that of germanium. A simple 
model for the number of counts in the full energy peak as a function of the 
y position of the beam, N(y), is given by: 



N(y) = N exp [- {fJ, aiPaiXai{v) + PCupCuX C u(y) + ^GePGe^Li (?/))] (1) 

where iVo is the number of full energy counts that would be observed in the 
detector without attenuation, px is the density of material X and xx{y) is 
the path length of the beam through material X and depends on the position 
of the beam, y. By fitting N(y) to data, x Li (y = 0) can be determined and 
identified with the minimum depth at which the full energy of the gamma 
rays is detected. This corresponds to the FCCD (full charge collection depth, 
defined above) which is traditionally identified as the dead layer thickness and 
determined by efficiency measurements. 

Fig. |3] shows data collected from a scan of the side of the LBNL PPC 
detector along the y direction fit to the function in equation [1] (solid line) 
to determine FCCD. At each position of the collimated source, the rate of 
events in the full energy photopeak at 59.5 keV was determined by integrating 
the energy spectrum near the peak (between 58keV and 61keV) and sub- 
tracting a background obtained from counts in the region above the peak, 
between 61 keV and 80keV to avoid any contribution from energy-degraded 
events in the peak. The unattenuated intensity, iVo, the position of the center 
of the crystal in relation to the collimator position, and FCCD were deter- 
mined from a fit performed using the MINUIT utilities in the ROOT software 



package from CERN [36|. Parameters, such as the attenuation coefficients, 



densities and positions of other materials (aluminum, copper) were fixed to 
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their known values. In order to account for a finite beam spot size, N(y) was 
convolved with a Gaussian. The fit to the data produced the best agreement 
when the standard deviation of the gaussian was equal to 2 mm. Systematic 
uncertainties from possible errors in parallax (to account for a possible mis- 
alignment of the beam with respect to the crystal coordinate system) were 
found to be the dominant systematic error in the determination of FCCD. 
To account for this, the angle between the y-axis of the collimator and the 
y-axis of the crystal was varied in the fit so that the uncertainty on FCCD 
from the fit includes this systematic contribution. The effect of changing the 
dimensions of the cryostat, the IR shield, and the germanium crystal was also 
examined but it was found that these dimensions are very well constrained 
by the fit and do not affect the uncertainty in determining FCCD. The full 
charge collection depth was determined to be 0.77 mm ± 0.12 mm. Since the 
dominant systematic parameters were also varied in the fit, the uncertainty 
includes systematic and statistical contributions. If the angle between the 
y-axes of the collimator and crystal was not allowed to vary, FCCD was de- 
termined to be 0.77 mm ± 0.06 mm, which is representative of the statistical 
uncertainty. The sharp decrease in the number of counts at y ~ 20 mm in 
Fig. [3] is due to the presence of a copper rod that is part of the crystal mount 
and was not incorporated into this simple analytical model. 

In order to verify the fit to the analytical function of equation [1] and 
confirm this model, a Monte Carlo simulation of the setup in Fig. [2] was per- 
formed using Geant4 version 4.9.2 37|. The geometry was reproduced in the 
model and a collimated beam of 59.5 keV gamma rays was simulated. The 
output of the Monte Carlo simulation was the distribution of energy deposi- 
tions (hits) in the crystal. Monte Carlo runs were performed by generating a 
fixed number of 59.5 keV gamma rays at each position of the collimator and 
source. For each position of the collimator, the distribution of the depth in 
the crystal at which events deposit their full energy was determined. This 
distribution is illustrated for two different y-positions (0 mm and -25 mm) 
of the collimator in Fig. HI Different depth profiles in the detector can be 
sampled by placing the source and collimators at different positions. In order 
to model the full charge collection depth, only events that deposit their full 
energy at a depth greater than FCCD (region 3 in Fig. H]) are considered to 
contribute to the full energy photopeak in the data. The number of events in 
the distribution illustrated in Fig. 0] above a given depth, FCCD, was used to 
build a probability distribution function (PDF) for the number of events in 
the photopeak as a function of the y-position of the collimator. The shapes 
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Figure 2: Projection of the setup used for the scan of the side of the detector, seen 
from below (z-axis into the page). The collimated beam of 59.5 keV gamma rays in the 
x direction is shown traversing the aluminum cryostat, the copper infrared shield, the 
lithium-diffused layer and the germanium. As the y-position of the collimator and source 
are varied, the path length of the beam through the different attenuating layers changes. 
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Figure 3: The rate of events in the 59.5 kcV photopcak as a function of the y-position of 
the collimated 241 Am source, from the setup illustrated in Fig. [5J The solid line shows a 
fit to the data using the the analytical function in equation [1] where the intensity of the 
beam, the center of the crystal, the thickness of the germanium dead layer, and the angle 
between the collimator and crystal y-axes were allowed to vary. The data were not fit 
beyond y = 20 mm, as the analytical model does not take into account the sharp fall off in 
events caused by a copper rod in the detector mount. The depth of the dead plus transition 
layer from the fit to the analytical function was determined to be 0.77 mm ± 0.12 mm. 
The dashed line shows a fit of the same data to the probability distribution functions 
obtained from a Monte Carlo simulation of the setup. The Monte Carlo model included 
the copper rod in the detector mount and reproduced the data beyond y = 20 mm. The 
depth determined using the Monte Carlo was 0.78 mm ± 0.08 mm, in excellent agreement 
with the analytical fit. 
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of the PDFs are very sensitive to the value of FCCD. These PDFs were built 
for a range of values of the full charge collection depth, as illustrated in Fig. 
E for FCCD= 0.2 mm, FCCD= 0.6 mm and FCCD= 1.2 mm. 




Figure 4: Normalized distributions of the depth at which Monte Carlo events have de- 
posited their full energy for two different y-positions of the collimated source. The col- 
limator position of y = mm corresponds to the center of the crystal and the position 
of y = —25mm is closer to the edge (Fig. [2]). The distributions show that by moving 
the collimated source, different depth profiles in the crystal are sampled. The area above 
0.78 mm (region 3) was determined by the fit in Fig. [3] and corresponds to events where 
the full energy is detected. The region between 0.48 mm and 0.78 mm (region 2) was 
determined from the fits in Fig. [10] and corresponds to a transition region where events 
have an energy degraded down to 26 keV and slow rise times. The region below 0.48 mm 
(region 1) corresponds to events where less than 26keV is detected and would correspond 
to a fully dead layer for a detector with a 26 keV threshold. 

A function that numerically interpolates between the PDFs was con- 
structed and used to fit the data in Fig. [3] using the MINUIT utility in 
the ROOT software package. The resulting fit is shown by the dashed line 
in Fig. [3] The PDFs for the different depths were normalized to have the 
same number of counts at y — and the fit function included an overall 
normalization to the number of events in the data as well as a parameter to 
allow the position of the center of the crystal to vary. Because the copper 
rod in the detector mount was included in the simulation, the range of the fit 
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Figure 5: Monte Carlo probability distribution functions (PDFs) of the number of full 
energy counts as a function of the collimator y-position for three different choices of FCCD. 
The PDFs were determined from the total number of counts above a given depth, FCCD, 
in the distributions shown in Fig. 0] (corresponding to region 3 in Fig. 2]). A function that 
interpolates between these PDFs was used to fit the data in Fig. [3l 
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was extended in the positive ^/-direction and reproduced the decrease in the 
number of counts from a copper rod in the detector mount at y ~ 20 mm. 

Systematic uncertainties were examined by varying the size of the col- 
limator in the Monte Carlo geometry as well as by tilting the beam below 
the x-axis and found to be negligible. Systematic effects from parallax were 
examined by letting the angle between the collimator y-axis and the detec- 
tor y-axis vary in the fit and were found to be the dominant contribution 
to the systematic uncertainty. With systematic and statistical uncertain- 
ties included, FCCD was determined to be 0.78 mm ± 0.08 mm, in excellent 
agreement with the fit using the analytical function. 



4. Evidence that slow, energy degraded pulses arise near the lithium- 
diffused surfaces 

Typical charge pulses from point contact detectors have a characteristic 
shape due to the sharply peaked weighting potential near the point contact. 
Pulses show an initial, slow, smaller amplitude component from the charge 
carriers drifting in a region of low weighting potential. This is followed by a 
larger, fast component from the short amount of time that the charge carriers 
spend in the region of high weighting potential near the point contact. A 
typical charge pulse from a 59.5 keV event in the PPC detector used in this 
work is shown in Fig. [6l along with its corresponding current pulse calculated 



with a 500ns averaging time using the method described in [29| (solid lines). 
Also shown in Fig. is a characteristically slow charge pulse along with its 
corresponding current pulse (dashed line). A simple metric to distinguish 
these types of events is either to measure the time that a charge waveform 
takes to rise from 10% to 80% or to use the asymmetry in the rise and fall 
times of the corresponding current pulse. 

The shape of the slow pulse cannot be explained by the drift of a single 
cloud of charge carriers to the point contact from any location in the crystal, 
as this would produce a charge pulse with a fast rise time. Such a long 
rise time indicates that the charge carriers arrive at the point contact over a 
period of time. 

The data collected by scanning the point contact surface with the colli- 
mated 241 Am source can be used to confirm the origin of the slow pulses. Fig. 
[7] shows a histogram of the distribution of the 10% to 80% charge waveform 
rise times as a function of the measured energy of the events in this scan. 



This is similar to Fig. 2 in Ref. |10J. The distribution shows a diagonal band 
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Figure 6: Comparison of the charge and current pulses for a 60keV (solid) and an energy- 
degraded 40 keV event (dashed). The current pulses were calculated from the charge pulses 
using the method described in [29I ] with an averaging of 500 ns. The lower energy event 
is a characteristic slow pulse that arises from the transition layer. It is easily selected 
based on either the rising time of the charge waveform or the asymmetry of the current 
pulse. The 10% to 80% rise time of the slow charge pulse is 1400 ns compared to 130 ns 
for the typical 60keV pulse illustrated here. For legibility in this figure, the waveforms 
were smoothed by a Gaussian filter. 
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of events originating from the 59.5 keV 241 Am peak with degraded energy and 
longer rise times. A box is drawn on the figure to highlight a clean sample 
of these events. 

Fig. [8] shows the fraction of events in this box relative to the number of 
events in the 59.5 keV photopeak as a function of the radial position of the 
collimated beam on the surface of the detector. The fraction of slow, energy 
degraded (SED) pulses rises dramatically near the edges of the crystal. This 
rise in the fraction of SED events occurs at the outer 5 mm, where the lithium- 
diffused wraparound is located (Fig. [1]). This result shows that significantly 
more SED pulses arise in the regions of a detector where a thick lithium- 
diffused contact is present. This is consistent with the interpretation in 



Ref. [10| that there exists a transition region in the lithium-diffused layer 



where partial charge collection occurs (degrading the energy), and that pulse 
shapes are slower as the charge carriers take time to diffuse into the fully 
active region of the crystal. 

Fig. [9]shows a comparison of the normalized 10% to 80% charge waveform 
rise time distributions of events with energies between 58 keV and 61 keV 
(near the 59.5 keV photopeak) (solid line) and that of events in the continuum 
between 20keV and 40keV (dashed line). There is an evident increase in 
the number of slow pulses at low energies, highlighting how events near the 
lithium-diffused layer can contaminate the lower part of the energy spectrum. 



5. Determination of the transition layer depth 

In Section [3j it was shown that different profiles for the depth of inter- 
actions in a detector could be sampled by scanning the side of the detector. 
The data from the scan of the side of the detector are used here to char- 
acterize the depth at which slow pulses arise. This section describes a fit 
of the fraction of slow pulses as a function of the collimator ^/-position to 
distributions obtained from the Monte Carlo simulation discussed in Section 

El 

It was found that, in the LBNL PPC detector, a slightly cleaner sample of 
the slow pulses could be obtained using the asymmetry of the current pulse 
(see Fig. E]) rather than the 10% to 80% charge pulse rise time. A cut to 
select slow pulses in the data was derived from the ratio of the 10% to 80% 
current rise time to the 80% to 10% current fall time. The fraction of slow 
pulses for each position of the collimator was defined as the number of slow 
pulses with energy between E m i n and 58 keV divided by the number of full 
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Figure 7: Distribution of the 10% to 80% rise times of charge waveforms as a function of 
the measured energy in the 241 Am source scan of the point contact face of the detector. 
A population of slow, energy degraded, (SED) pulses is evident as the diagonal band 
emerging from the 59.5 keV photopeak. A box is shown to select a clean sample of these 
events; most of the data (83%) lie in the band below 300 ns. 
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Figure 8: The fraction of SED pulses in the scan of the point contact surface of the detector 
as a function of the radial position of the collimated beam (the beam points upwards when 
referring to Fig. Q] and a radius of zero corresponds to the point contact). For a given 
position of the collimator, the fraction of slow pulses is taken as the number of events in 
the box shown in Fig. [7] divided by the number of events in the full energy photopeak 
at 59.5 kcV. The fraction of SED pulses increases dramatically at the edge of the crystal, 
near a radius of 25 mm, where the 5 mm-wide lithium wraparound is located, confirming 
that these are related to charge depositions near the lithium-diffused region of germanium. 
The diameter of the beam spot was approximately 3 mm in these data. 
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Figure 9: Normalized distribution of the 10% to 80% rise times of the charge waveforms 
for events with energies in the 59.5 keV photopeak (solid) and within the continuum be- 
tween 20keV and 40keV (dashed). Interactions in the lithium-diffused layer result in a 
population of characteristically slow pulses contaminating the lower part of the energy 
spectrum. 
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energy events in the 59.5 keV photopeak (determined in the same way as for 
the measurement of FCCD). The number of slow pulses between E min and 
58keV was corrected by subtracting an energy-independent background of 
slow pulses estimated by integrating the number of slow pulses in a range from 
61 keV to 80 keV. It was verified that the inferred statistical uncertainty from 
the background subtraction was consistent with data taken with no source 
present. The fraction of slow pulses as a function of the collimator position 
is shown in Fig. [10] for E min =18keV and i? m j n =38keV. As the the threshold 
for the minimum energy is lowered, the fraction of slow pulses increases, as 
expected from Fig [3 
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Figure 10: The fraction of SED events in the scan of the side of the detector as a function 
of the y-position of the collimator (in the same geometry as in Fig. [2]). The fraction 
of slow pulses in the data is taken as the number of slow events with measured energies 
between E min and 58keV divided by the number of events in the 59.5 keV photopeak. 
The distribution of slow events is shown for values of £' m i n =18keV and E m i n =38keV. 
The data with the lower energy threshold have a larger fraction of SED events. These 
are fit to PDFs obtained from the Monte Carlo simulation (solid lines) to determine the 
depth at which slow pulses with an energy above E m i n are detected. It is worth noting 
that the Monte Carlo PDFs did not need to be normalized to the data and that the depth 
of the transition layer can be determined at any position from the ratio of slow pulses to 
full energy events independently from the source strength. 

The distributions of slow pulses as a function of the collimator y-position 
for different values of E m i n in the data were then fit to a set of PDFs derived 
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from the Monte Carlo. The PDFs were obtained in a way similar to those 
that were used to determine the full charge collection depth. In this case, 
for each position of the collimator in the Monte Carlo data, the fraction 
of slow pulses is taken as the number of events that deposit their energy 
between the transition layer depth, TLD, and FCCD divided by the number 
of events that deposit their full energy at a depth greater than FCCD (see 
Fig. HJ). TLD(E min ) is identified as the depth at which full energy depositions 
are detected with an energy above E m i n . The denominator is simply the 
number of full energy events that was used previously to determine FCCD. 
The fraction of slow pulses in the Monte Carlo data corresponds to the ratio 
of the area in region 2 to the area in region 3 in Fig. SJ 

A set of Monte Carlo PDFs were built with different values of the TLD. 
A function was then constructed to interpolate between these PDFs and fit 
the data in Fig. [10] to determine the value of the transition layer depth. 
The position of the crystal and uncertainty from parallax were determined 
from the fits for the full charge collection depth in Fig. [3] and allowed to 
vary within one standard deviation of those results. The fits were performed 
using different values of E min in the data. The solid lines in Fig. [TU] show 
the results of this fit to the data with E , min =18keV and i? m i n =38keV. Data 
were only fit for negative values of the ^/-position of the collimator to remove 
any effect that may be due to the copper rods in the mount. It was verified 
that the determined value of TLD was not affected when the fit range was 
extended, although the quality of the fit worsened. 

The dominant systematic uncertainty in determining TLD was the full 
charge collection depth. The two quantities are directly correlated as can be 
seen in Fig. [Til which shows the value of TLD that was determined for data 
with E min =2QkeV as a function of FCCD. The statistical and systematic 
errors from the fit are are too small to see in Fig. [TT] and are comparatively 
negligible. The uncertainty on TLD was determined by using the range 
corresponding to FCCD = 0.78 mm ± 0.08 mm, as determined in Section [3] 
and is illustrated by the shaded region in Fig. [TTJ The depth at which slow 
pulses with a minimum detected energy of 26 keV arise was determined to be 
0.48 mm ± 0.08 mm. The location where these slow pulses arise corresponds 
to region 2 in Fig. [H 

The analysis was performed for different values of E m i n ranging from 
18keV (the energy threshold that was used to collect the data) to 54keV in 
4keV increments. The value of TLD that was determined as a function of 
the minimum energy of the slow pulses is shown in Fig. H2J The detector, 
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Figure 11: Transition layer depth (TLD) for slow pulses with £" m i„=26keV determined 
from a fit similar to those in Fig. QJjas a function of the full charge collection depth (FCCD) 
that was used to generate the PDFs. The horizontal shaded area shows the uncertainty 
in TLD that is determined from the allowed range in FCCD obtained from the fit in Fig. 
(shown with the vertical lines). This is the dominant systematic uncertainty in the 
transition layer depth. FCCD was determined to be 0.78 mm ± 0.08 mm in Section [3] and 
the resulting value for TLD is thus 0.48 mm ± 0.08 mm, for slow events with an energy 
above 26 kcV. 
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although capable of a much lower energy threshold, was only operated down 
to 18keV because of an instability in the preamplifier and the trend is not 
demonstrated to lower energies. Nonetheless, it is clear that high energy 
events that deposit their energy in the transition region result in slow pulses 
with a degraded measured energy. As the depth of an energy deposition 
increases, more energy is detected. Fig. 2 in Ref. 10] as well as data taken 
by the the Majorana collaboration with the detector described in Ref. 23 
show that the slow, energy-degraded, events from a 59.5 keV 241 Am gamma 
ray source extend to below 1 keV. This is of significant concern to experiments 
operating with a very low energy threshold. The shape of the curve in Fig. 
12 at low energies is unknown, and is likely different for different detectors, 
depending on how the n+ contact is fabricated 



Q_ 
CD 
"D 
i 

CD 

>, 

CO 



c 
cfl 




10 



20 30 40 50 60 

Min. energy of slow pulses (E ) [keV] 



Figure 12: Transition layer depth as the minimum energy of the SED pulses is varied. 
As lower energy SED events are considered, the transition layer depth decreases. The 
truly dead layer of the detector, where no energy depositions would be detected, may be 
quite small, as the data are extrapolated to lower energies. Understanding this curve for a 
detector is critical in estimating the background from slow pulses that one would observe 
at a very low energy threshold. 



In order to understand this result, a similar set of analyses were per- 
formed using the detector described in Ref. 23|. This is a modified "Broad 
Energy Germanium" (BEGe) detector fabricated by CANBERRA and is 
similar to the one used by the CoGeNT collaboration. The modifications 
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include a thicker lithium-diffused layer, a smaller point contact (4 mm di- 
ameter) compared to the commercially available BEGe detectors, and a low 
background copper cryostat. The detector is operated as part of the Ma- 
jorana Low-background BEGe at the Kimbalton Underground Research 
Facility (MALBEK) experiment. 

The full charge collection depth for the MALBEK detector was obtained 
with a method similar to that described in Ref. jsj using a 133 Ba source. 
An energy spectrum was collected using a 0.9 /iCi 133 Ba calibration source 
placed 25 cm from the end-cap of the detector (side opposite of the point 
contact). The ratio of the number of events in the 81keV photopeak to that 
in the 356 keV photopeak was then determined. The experimental setup and 
cryostat were simulated using the Geant4-based MaGe simulation software 



381 ] . A full charge collection depth was implemented in the post-processing of 
the simulation and was adjusted so that the ratio of the number of counts in 
the two photopeaks matched the one measured in the data. Fig. [TBI shows the 
ratio from the simulation as a function of FCCD (points with statistical error 
bars) as well as the ratio obtained from the data (horizontal lines). FCCD 
for this detector (vertical lines in Fig. [T3"j) was determined to be 0.933 mm 
± 0.018mm (stat.) ± 0.092 (sys.). In determining the systematic error, 
uncertainties in the length of the crystal and the thickness of the copper 
end-cap in the cryostat dominated over other effects such as uncertainties in 
the relative intensities of the photopeaks, live time, and the determination 
of the number of counts in the peaks. 

In order to probe the transition layer in the MALBEK detector, the 
end-cap of the detector was illuminated with an un-collimated 10 /LtCi 241 Am 
source of 59.5 keV gamma rays. The experimental setup was simulated and 
the transition layer depth was modeled in the same way as for the LBNL 
detector. SED pulses were selected based on the 10%-90% rising time of the 
charge waveform and backgrounds were determined using data taken with 
no source present. Fig. [TH shows the transition layer depth determined for 
the MALBEK detector as the minimum energy of the SED pulses is varied 
down to 5 keV. The dominant systematic uncertainty in determining TLD for 
this detector was also found to be the uncertainty in determining FCCD (see 
above), which resulted in an 0.094mm systematic uncertainty included in 
the error bars in Fig. [TH The shape of the curve in Fig. [THis similar to that 
in Fig. [12] and does not show any discontinuities over this energy range. The 
derivative of the two curves are different and suggest that the widths of the 
transition region and the fully dead layer are different for the two detectors 
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Figure 13: Determination of the full charge collection depth for the MALBEK detector 
using data from a 133 Ba calibration source. The points with statistical error bars show the 
ratio of the number of counts in the 81keV and 356 keV photopcaks from a simulation of 
the detector as a function of FCCD. The horizontal lines show the ratio that was measured 
in the data and the vertical lines show the value of FCCD determined by comparing the 
data and the simulation. For this detector, FCCD was determined to be 0.933 mm ± 
0.018mm (stat.) ± 0.092 (sys.). 
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and likely depend on the process used to fabricate the lithium- diffused layer. 
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Figure 14: Transition layer depth as the minimum energy of the SED pulses is varied in 
the MALBEK detector. Data were recorded down to 5 keV and show that this curve has 
a smooth behavior down to these low energies. The dominant systematic error shown on 
these data points is from the uncertainty in determining FCCD. 



The results presented in this work show that one can characterize the 
depth of the transition region and that a straightforward Monte Carlo simu- 
lation of the depth of energy depositions can be used to estimate the level of 
contamination from slow pulses. Although a very detailed scan of the LBNL 
PPC detector was performed, the transition layer depth and the full charge 
collection depth can be determined with few measurements. For example, 
FCCD can be determined using standard photopeak efficiency measurements. 
A particularly useful technique is to use gamma rays emitted from a 133 Ba 
source and to compare the ratio of events in the low energy phototopeaks 
(around 81keV) to the number of events in the dominant line at 356 keV to 
obtain FCCD, as was done for the MALBEK detector. These data could 
also be used to determine TLD by measuring the ratio of the number of slow 
pulses that come from the lines at 81 keV to the number of full energy events 
in the corresponding photopeaks. A flood source measurement with 241 Am 
was used for MALBEK to determine the transition layer depth on the surface 
opposite to the point contact. It is worth noting that these measurements 
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only require a detailed knowledge of the absorbers between the calibration 
source and the detector and that they are independent of the calibration 
source strength. 

The level of detail required to fully characterize a detector and quanti- 
tatively predict the number of SED pulses depends on an understanding of 
the uniformity of the lithium-diffusion. For the LBNL PPC in this work, the 
measured distributions of full energy events (Fig. and of SED pulses (Fig. 
[TO]) could be fit over the entire range of ^-positions of the collimator sug- 
gesting that the lithium-diffused layer along that axis is uniform. Although 
FCCD and TLD can be determined by a measurement at a single position 
on the detector, experiments requiring a thorough understanding of the SED 
contributions to their spectrum should perform a measurement to quantify 
the uniformity of the lithium-diffused contact. This can be achieved with a 
combination of flood source measurements, as was done for the MALBEK 
detector, and measurements performed using a collimated source. 

6. Conclusion 

This paper presented data taken with the MALBEK detector as well 
as with a custom PPC detector fabricated at Lawrence Berkeley National 
Laboratory (LBNL). The mount holding the LBNL crystal and the geometry 
of the n+ contact of this system allowed effects from low energy gamma rays 
interacting in the lithium-diffused outer n+ contact of the detector to be 
examined in detail. 

The existence of a transition layer between the active region of the de- 
tector and the outer n+ lithium-diffused contact was confirmed. Charge 
depositions in this transition layer were observed to result in partial charge 
collection (energy degradation) as well as characteristically slow rise times 
of the corresponding charge pulses. A method was introduced to directly 
measure the full charge collection depth as well as the depth at which the 
transition layer begins. It was shown that energy depositions near the sur- 
face of the detector contacts can still result in detectable amounts of energy 
and that the truly dead layer can be smaller than is typically assumed. Data 
taken with the MALBEK detector suggested its transition and fully dead 
layers had different widths than those of the LBNL detector, indicating that 
the process for fabricating the lithium- diffused contact may play an impor- 
tant role. Detailed spectral simulations of germanium detectors with thick 
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lithium- diffused contacts need to account for these effects in order to repro- 
duce spectral features at energies below known photopeaks. 

The slow, energy-degraded, (SED) pulses are of particular concern to ex- 
periments using germanium detectors with a very low energy threshold, such 
as for dark matter detection. Matters are complicated by the fact that cuts 
to remove these events based on pulse shapes (such as those used in Ref. [lol ] ) 
become difficult to implement at energies of a few keV, where pulses are noisy 
and calculations must rely on effective de-noising techniques. Experiments 
using data from the lowest energies must demonstrate that contributions 
from SED events are understood. This is particularly important if there are 
gamma or x-ray lines whose SED events contribute to the continuum and 
a detailed knowledge of the spectrum is required. This work presented a 
method to characterize energy depositions in the lithium-diffused contacts 
which can be used to quantify the contributions to backgrounds from SED 
pulses. 

Mitigation strategies for experiments looking for signals at lower ener- 
gies should include detailed characterizations of the detectors and of cuts 
based on pulse shapes, along with the demonstration that the associated 
systematic uncertainties are well understood. Detector development work 
and research into fabricating uniform and/or thinner hole-blocking contacts 
on the detectors, which could result in a sharp transition layer, are another 
research avenue. Promising work has been done to show that one can make 
thin, effective, hole-blocking contacts using amorphous silicon or germanium 



31| and yttrium |39|. Furthermore, n-type point contact detectors [1[ may 
have similar performance as p-type detectors at low energies; these can be a 
candidate for low energy experiments as they would almost entirely remove 
this source of background. However, if detectors with thin contacts are used, 
the benefit of shielding low energy backgrounds with the thick traditional 
contact of p-type detectors would be lost. 
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